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1. Introduction 

Tetrahedrally coordinated binary semiconductors of A^B^~^ type having total of 8 
valence electrons are compounds of great technological importances. They have been 
widely studied both theoretically and experiment ally [H E]. An interesting aspect of 
these compounds is that they can belong to either cubic or in hexagonal symmetry, 
namely in zincblende(ZB) or in wurtzite(WZ) structure. The difference between the 
two structures is subtle, both being of tetrahedral coordination, and also the associated 
total energy difference is small, of the order of few tens of meV/atom. These particular 
aspects have made the ZB vs WZ structural stability issue a topic of significant interest 
that have been discussed in literature both in terms of model calculations as well as 
in terms of first principle calculations P lll[5l|6l[7l|8l|9l[Tni[IIl[l2]. Phillips and Van 
Vechten [5] first gave a well defined prescription in terms of ionicity to predict the 
crystalline phase in binary octet semiconductors of A'^B^~'^ type, with the help of 
the dielectric theory. While their theory nicely described the transition from four-fold 
to six-fold coordination and possibly also to eight-fold coordination, it was not able 
to distinguish clearly between structures having the same coordination. P. Lawaetz[7] 
investigated the stability of WZ structure more closely. Taking into account the short- 
range elastic forces and long-range Coulomb forces, he correlated the ZB vs WZ stability 
with the deviation of the axial ratio - from the ideal value, which was further related 
with a structure-independent charge parameter. However, the results were not fully 
satisfactory because of limited knowledge about the long-range Coulomb-effects in 
partially ionic materials. Christensen et al. [H] studied the structural phase stability 
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of 34 elemental and compound semiconductors from first principle electronic structure 
calculations which involved mostly four-fold to six-fold coordination transitions. They 
studied the chemical trends by calculating the valence charge densities using linearized 
muffin-tin orbital (LMTO) based localized basis sets. Their calculated ionicities from 
tight-binding representation of LMTO Hamiltonian, gives better chemical trend than 
Phillips value. Recently, the ZB - WZ polytypism in binary octet semiconductors have 
been studied by the group of Alex Zunger[T2] in terms of quantum mechanically defined 
atomic-orbital radii. They showed a linear scaling between the ZB - WZ energy difference 
and renormalized orbital radii and successfully predicted the exact structural trends in 
most of the octet compounds, except some few cases. The general understanding that 
emerged out of these calculations is that it is the competition between the covalency 
and the ionicity effects that determines the relative stability of ZB vs WZ structures, 
with covalency favoring ZB structure and ionicity favoring wurtzite structure. 

Experimentally, the stable crystal phase of CdS is hexagonal wurtzite structure [13]. For 
CdSe, zincblende is the stable low temperature phase and above a critical temperature, 
it transforms to wurtzite structure [H]. On the other hand, CdTe always stabilizes in 
cubic zincblende structure [15]. However, there is a few indication of metastable wurtzite 
growth for CdTe in some special situation [16]. In the present work, we revised the 
issue of relative stability of ZB and WZ in the CdX series, with X = S, Se, Te using 
Nth-order muffin-tin orbital (NMTO) technique|17j within the framework of density 
functional theory (DFT). In particular, we employed the "direct generation of Wannier- 
like orbitals" feature of NMTO technique for attaining the microscopic understanding 
in this context. 

2. Crystal structure 

The ZB structure consists of two interpenetrating face centered cubic (FCC) sublattices, 
one of atom A, the other of atom B, displaced from each other along the body diagonal 
by |, a being the lattice constant for the ZB structure. On the other hand, an ideal WZ 
structure consists of two interpenetrating hexagonal closed packed(IICP) sublattices, 
one of atom A, the other of atom B, displaced from each other by |c along the c-axis. 
These result into two different stacking sequences: ABCABC. . . along [111] direction for 
ZB and ABAB. . . along c-axis for WZ. The ZB structure corresponds to the staggered 
conformation of atomic arrangement along [111] body diagonal, while the wurtzite 
structure corresponds to the eclipsed conformation as seen looking down the c-axis. 
The nearest neighbor (tetrahedral bond) arrangements in the ZB structure and in the 
ideal WZ structure are identical. The main difference starts to come in the relative 
position of 3rd nearest neighbors and beyond. Also the arrangement of the distant 
atoms along the four different tetrahedral bonds are different for a WZ structure. The 
structural differences between these two phases have been shown in Fig. [1] 

The lower symmetry of the WZ structure allows for a distortion with the c/a ratio 
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wurtzite zincblende 




Figure 1. (a) and (b) show the crystal structure in WZ and ZB symmetries. Two 
different colored atoms denote cation and anion, (c) WZ structure shows the ABAB. . . 
stacking of atoms, while (d) ZB structure shows the ABCABC. . . stacking. Eclipsed 
and staggered conformation of atoms along c-axis in WZ structure and along [111] body 
diagonal in ZB structure have been shown in (e) and (f) respectively. 

deviating from the ideal value of = 1.633. A further kind of distortion, that 
is possible in the hexagonal structure and forbidden in the cubic, is that the two 
interpenetrating sublattices are slightly displaced from each other so that an atom of one 
sublattice is no longer at the centroid of the other sublattice. In order to capture these 
distortions, it is necessary to define another internal parameter "u" , which is the nearest 
neighbor distance along c-axis, in unit of c. The value u in the ideal WZ structure is 
I = 0.375. Distortion of u from its ideal value gives rise to a dipole at the center of 
each tetrahedron, containing an atom of opposite kind, the magnitude of which depends 
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upon Am [Am = u - 0.375] and the ionicity of the chemical bonds. This effect also lowers 
the ionization energy. The structural parameters of CdS, CdSe and CdTe used in the 
present calculation are listed in Table [TJ 

Table 1. Structural parameters of ZB and WZ phases for CdX (X = S, Se, Te) series. 
Lattice parameters have been taken from Ref. [12] for both ZB and WZ phases of CdS, 
from Ref. [18] for ZB CdSe, from Ref. [19] for WZ phase of CdSe, from Ref. [15] for 
ZB CdTe and from Ref. [16] for WZ CdTe. 



Compounds 






WZ 






ZB 




a(A) 


C 

a 


u 


A(f) 


Am 


a(l) 


CdS 


4.121 


1.621 


0.377 


-0.012 


0.002 


5.811 


CdSe 


4.299 


1.631 


0.376 


-0.002 


0.001 


6.077 


CdTe 


4.57 


1.637 


0.375 


0.004 


0.000 


6.492 



It shows that the deviation A(f ) [A ) = f - 1-633] is -ve and largest for CdS, while it 
is -|-ve for CdTe. The deviation from the ideal c/a ratio is tiny for CdSe. The deviation 
Am is found to be opposite in sign to that of A(^). A part of the deviation Am is 
expected to come due to nearest neighbor bond bending or bond stretching caused by 
A(^). The limits on this part, as calculated by Keffer and Portisjl] for AIN, are u 
= 0.380 to maintain equal nearest neighbour bond lengths and u = 0.372 to maintain 
equal four tetrahedral angles. However, measured value of m for AIN is 0.385[1]. This 
means there must be additional contribution in the distortion of u from ideal value, 
which is coming from the long-range Coulomb interaction. NMTO technique based on 
full self-consistent DFT calculations includes both the effect of short-range as well as 
long-range interaction in providing a complete understanding. 

3. Total energy calculations 

In order to check the relative stability of the ZB and WZ structures of the CdS, 
CdSe and CdTe, we carried out total energy calculations using the energetically 
accurate pseudo-potential basis set. The calculations have been performed with ultra- 
soft pseudopotentialpO] and local density approximation for the exchange-correlation 
functionals implemented in the Vienna Ab-initio Simulation Package(VASP) [2T| . The 
wavefunctions are expanded in the plane wave basis set with a kinetic energy cut- 
off of 250 eV which gives convergence of total energy sufficient to discuss the relative 
stability of various phases. The calculations have been carried out with a k-space grid 
of 11 X 11 X xll. The obtained total energies and the band gaps are listed in Table 
|2l As found, the total energy differences are indeed very tiny; of the order of few meV. 
Our calculations correctly show that CdS to be stable in WZ phase and CdTe to be 
stable in ZB phase. CdSe which is considered as a boarder line case is found to be 
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stable in ZB phase in agreement with the pubhshed results [T^. The band gaps in WZ 
and ZB structures, both being direct band gaps for CdS, CdSe and CdTe are similar for 
a given compound. It is found that the band gap decreases monotonically in moving 
from CdS to CdTe in both ZB and WZ symmetries. This monotonic decrease in band 
gap indicates the enhanced metallicity across the CdS-CdSe-CdTe series which points 
to CdTe being more covalent compared to CdS. 

Table 2. Total energy per formula, energy difference AE — E{ZB) - EiWZ) and 
band gaps. + ve sign of AE indicates WZ is more stable than ZB and vice versa. 
The numbers within the parenthesis in the columns of band gap, correspond to their 
experimental values taken from Ref. (35] 



Compounds 


Energy (eV/formula) 


AE (mcV) 


Band ; 


?ap(eV) 




ZB 


WZ 




ZB 


WZ 


CdS 


-6.8650 


-6.8665 


1.5 


0.94(2.55) 


1.02(2.58) 


CdSe 


-6.7548 


-6.7523 


-2.5 


0.77(1.90) 


0.72(1.83) 


CdTe 


-6.0248 


-6.0221 


-2.7 


0.69(1.60) 


0.65(1.60) 



4. Calculation of ionicity 

As already pointed out, the relative stability between ZB and WZ phases of B^~^ 
semiconductors is dictated by the competition between the covalent bonding and the 
electrostatic energy given by the difference between the cation and anion energy levels. 
In a homopolar semiconductors like Si, it is only the first term that survives. For 
a tetrahedrally coordinated semiconductor, it is most natural to think in terms of sp^ 
hybrids. Considering the sp^ hybrid energy as EgpZ = ^"'^^^p the energy level separation 
between cation and anion is given by AEgpS = E'^^s — E'^^s ; where c and a denote cation 
and anion respectively. The hybridization contribution on the other hand is related 
to the hopping integral between the cation and anion sp^ hybrids. It is given by h 
= J l^s"- + Px + Py + Pt\H\s'^ — pI — Py — pI^ considering the [111] bond of ZB structure 
and h = ^ {s"" + p'^\H\s'^ ~ pl) considering the [001] bond of WZ structure where H is 
the tight binding Hamiltonian in the sp basis of cation and anion. The ionicity is then 
defined as 

(AE,,3)2 + (2/z)2 ^ ^ 

In order to extract the energy level separation and the hopping interaction, we employed 
the NMTO-downfolding technique. In this method, a basis set of localized orbitals is 
constructed from the exact scattering solutions (partial waves and Hankel functions) 
for a superposition of short-ranged spherically-symmetric potential wells, a so called 
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muffin-tin (MT) approximation to the potential. The basis set is constructed from 
the scattering solutions at a mesh of energies, eo, ei, cat. At those energies, the 
set provides the exact solutions, while at other energies, the error is proportional 
to (£■ — eo)(-E' — ci) . . . (i? — eAf). The basis set of NMTOs is therefore selective in 
energy. Moreover, each NMTO satisfies a specific boundary condition which gives it 
a specific orbital character and makes it localized. The NMTO basis set due to its 
energy selective character is flexible and can be chosen as minimal via the downfolding 
procedure which integrates out the degrees of freedom that are not of interest. The 
downfolded NMTO set spans only the selected set of bands with as few basis functions 
as there are bands. For the isolated set of selected bands the NMTO set spans the Hilbert 
space of the Waannier functions, that is, the orthonormalized NMTOs are the Wannier 
functions. The Wannier orbitals are therefore generated directly in this method, which 
may be contrasted with the techniques where Wannier functions are generated out of the 
calculated Bloch functions as a post processing step [23]. The downfolded Hamiltonian 
in the Wannier basis provides the estimates of onsite energies and the hopping integrals. 

In order to compute the ionicity parameter defined in Eqn. (1), the NMTO calculations 
have been carried out with Cd(spc?) and X(sp) basis. NMTO calculations being not yet 
implemented in self-consistent form, the self consistent calculations have been carried 
out in the LMTO method. The self consistent MT potential out of these LMTO 
calculations has been used for the constructions of NMTOs, which in the present case is 
the standard LMTO all-electron CdXEE' atomic spheres potential for the ZB structure 
and Cd2X2E2E2 for the WZ structure. E and E' are two different empty spheres used 
to fill the space. The calculated estimates of AEgpS, hopping integral h and ionicity 
parameter /j for CdX series are listed in Table [3l The NMTO approach successfully 
brings out the right trend within the CdX series both for the ZB and the WZ structures, 
i.e CdS has the maximum ionicity and CdTe has the least. For comparison, the results of 
previous calculations by Christensen et al. |Tl] using the tight-binding LMTO approach 
and that by Phillips [2] using dielectric theory are also listed. As it is found, while 
NMTO could bring out the right trend within the fine differences, the two previous 
approaches could not capture it properly. 



Table 3. Covalent gap AEgpS, hopping term Fih ~ -2h and ionicity fi for the CdX 
series in ZB and WZ structure. fi{C) means calculated value of ionicity by Christensen 
et al. [n] and /^(P) that of Phillips taken from Ref.;2j 



Compounds 


AE,p3 


(eV) 


Eh = 


-2h (eV) 


/. 




/.(C) 


MP) 




ZB 


WZ 


ZB 


WZ 


ZB 


WZ 


ZB WZ 


ZB WZ 


CdS 


7.13 


7.17 


6.37 


6.29 


0.556 


0.565 


0.794 . . . 


0.685 . . . 


CdSe 


6.70 


6.69 


6.16 


6.04 


0.542 


0.551 


0.841 . . . 


0.699 . . . 


CdTe 


5.36 


5.45 


5.79 


5.74 


0.462 


0.475 


0.739 . . . 


0.717 ... 
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5. Microscopic understanding in terms of calculated Wannier functions 

In order to gain further insights, we went a step ahead and constructed the truly minimal 
NMTO sets with the sp-orbitals placed exclusively on the anion(X) site and downfolded 
all the orbitals at the cation(Cd) site, except the cation d orbitals. This gives rise to the 
basis with only 4 sp-orbitals out of 8 sp-orbitals in the ZB unit cell and 8 sp-orbitals out 
of 16 sp-orbitals in the WZ unit cell. The energy points are chosen in the way so as to 
span only the valence bands. The comparison of the downfolded valence-only bands and 
the full band structure can be made as good as possible by making the energy mesh finer 
and finer. The plot shown in Figj2]for CdS with the choice of 4 energy points already 
shows the downfolded bands to be indistinguishable from the full band structure in the 
scale of the plot. Similar agreements are found also for CdSe and CdTe. 



Cd(d)S(sp) - Zincblende Cd(d)S(sp) - Wurtzite 




L G XWLK G G KUGALH A 



Figure 2. Downfolded bands (thick line) with the basis sets where sp-MTOs are on 
S atom and d-MTOs are on Cd atom, are compared with NMTO all bands (thin line) 
both in ZB(left) and WZ(right) structure of CdS. The energy mesh for downfolding 
have also been shown. 

FiglH] shows the plot of the orthonormalized p-NMTO (N=3) centered at the anion site 
and pointing to the neighboring Cd site along the [111] direction for the ZB structure 
and [001] direction for the WZ structure. As it is clearly seen, the red lobe at the left- 
hand side which is a mark of the covalency effect, systematically increases in moving 
from S to Se to Te. Had it been plotted for the case of homopolar system like Si which 
is a perfectly covalent compound, the plot would have been perfectly symmetric with 
the red lobes being symmetric between left hand and the right hand side. These plots 
reconfirm the conclusion that the ionicity decreases and the covalency increases across 
the CdX series. 

Fig. m shows the plots of directed p-NMTO same as in Fig. [3] but for four different 
tetrahedral bond directions of WZ and ZB structure for CdS. While the four bond- 
centered p-NMTOs look identical for ZB structure, the bond-centered p-NMTO directed 
along the [001] direction of the WZ structure, looks different from the rest of the three. 
The difference primarily comes from the tail sitting at Cd position at c{l—u) measured 
from the central anion position at X(0, 0, 0) along the [001] axis ( see Fig. [5]). There is 
no equivalent neighboring atom along three other directions. 
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Figure 3. Tlic contour plots sliow the anion p[lll] MTO in ZB structure (left panels) 
and anion p[001] MTO in WZ structure (right panels) of CdX (X = S, Se, Te). The 
top most panels correspond to CdS and the bottom most to CdTe. In each case, 35 
contours have been drawn in the range -0.15 to 0.15 electrons/Bohr'^. From top to 
bottom, the ionicity decreases and the covalency increases. 



In order to understand this effect, we have carried out model calculations where ^ j 
ratios have been made -0.1, -0.05 and 0.1. The corresponding "u" parameter in each 
case has been obtained by total energy minimization fixing the lattice parameters. The 
deviation in "u", Au is found to roughly obey the relationship with A(^) as Am = - 

(lis) ^ (a) where ^, the bond-bending parameter, is 2.0, given by P. Lawaetz[7]. We 
computed the ionicity parameter for each of the model systems by NMTO downfolding 
technique. In Table H] we list the ionicity parameter for the model WZ systems, together 
with actual CdTe case. As it is evident the -ve deviation from ideal (^^^ ratio makes the 

compound more ionic as pointed out previously. Changing A (^^^ ratio from -ve to -|-ve 

passing through the ideal A (^^^ ratio, the ionicity decreases and covalency increases. 
This is in agreement with the structural parameters of CdS, CdSe and CdTe. CdS being 
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the most ionic among the three, shows the largest negative deviation in (^^^ . (^^^ ratio 
for CdSe is close to ideal, while deviation in (^^j ratio for CdTe becomes positive. 

Table 4. Covalent gap AEgpS, hopping term ~ -2h and ionicity fi for model CdTe 
in WZ structure, with A and u as given in the first column. 





CdTe 




AE,p3 (eV) 


Eh = -2h (eV) 




A(f) 


= -0.100 and u = 


0.4044 


5.49 


5.33 


0.515 


A(f) 


= -0.050 and u = 


0.3903 


5.48 


5.57 


0.491 


A(f) 


= +0.004 and u = 


= 0.375 


5.45 


5.73 


0.475 


A(f) 


= -HO. 100 and u = 


: 0.3432 


5.45 


6.13 


0.441 



In Fig. ini we show the downfolded p-NMTO for the model WZ systems together with 
actual WZ CdTe system along [001] direction and one of the three other directions, 
namely [-1, 0, ]. As it is evident, changing ^ (f) from -ve to -|-ve makes the 
two p-NMTOs directed along the vertical [001] bond and one of the other directions, 
looking alike, diminishing the tail effect siting at Cd[0, 0, c(l - u)] site. This is driven 
by the covalency effect which prefers an isotropic arrangement. However, it cannot be 
completely achieved within a hexagonal symmetry, a compound with nearly ideal (^^j 

or positive A (^^^ ratio therefore prefers to stabilize in ZB symmetry, satisfying the 
isotropic distribution completely. 

6. Summery 

Using NMTO-downfolding technique we have revisited the problem of ZB vs WZ 
symmetry in case of A^B^''^ semiconductors. In particular, we have considered the 
CdX series with X =S, Se, Te. Our computed ionicity factors using accurate NMTO- 
downfolding successfully brings out the right trend within the CdX series - CdS being 
most ionic stabilizes in WZ symmetry while CdSe and CdTe being more covalent 
stabilizes in ZB symmetry. Our NMTO constructed Wannier functions corresponding 
to only valence bands provide nice demonstration of this fact. The tendency towards ZB 
stability is governed by the covalency which prefers isotropic nature of the tetrahedral 
bonds. 
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Figure 4. The contour plots show the p-MTO of S along the four nearest neighbor 
tetrahedral directions in ZB structure (left panels; from top to bottom they are p[lll], 
p[lli], p[ni]) and in WZ structure (right panels; from top to bottom they are 



p[001],p[l,0, ^f],p 
are same as in Fig. [3] 
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Figure 5. (a) Positions of atoms along [001] direction of wurtzite structure, (b) 
positions of atoms in unit of lattice constant along [111] direction of zincblende 
structure. 
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Model CdTe in wurtzite structure y j. 
Te - p[0, 0, 1] Te - p[ 1, 0,^ — ] 




Figure 6. The above contour plots show the p-MTO of Te in WZ structure of model 
CdTe. The top panels correspond to p[001] and bottom panels correspond to p[l, 0, 
-^^]. From left to right ionicity decreases and covalency increases and distribution 
becomes more isotropic. The contours chosen are same as in FiglS] 



